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ABSTRACT 

We present deep Spitzer mid-infrared spectroscopy, along with 16, 24, 70, and 850 /im photometry, 
for 22 galaxies located in the Great Observatories Origins Deep Survey-North (GOODS-N) field. The 
sample spans a redshift range of 0.6 ^ z < 2.6, 24 /im flux densities between ~0.2— 1.2 mJy, and 
consists of submillimeter galaxies (SMGs), X-ray or optically selected active galactic nuclei (AGN), 
and optically faint {zab > 25 mag) sources. We find that infrared (IR; 8 — 1000 /im) luminosities 
derived by fitting local spectral energy distributions (SEDs) with 24 /xm photometry alone are well 
matched to those when additional mid-infrared spectroscopic and longer wavelength photometric 
data is used for galaxies having z < 1.4 and 24 /j,m-derived IR luminosities typically ;< 3 x 10^^ Lq. 
However, for galaxies in the redshift range between 1.4 ^ z ^ 2.6, typically having 24 ^m-derived 
IR luminosities > 3 x 10^^ Lq, IR luminosities are overestimated by an average factor of '^5 when 
SED fitting with 24 fim photometry alone. This result arises partly due to the fact that high redshift 
galaxies exhibit aromatic feature equivalent widths that are large compared to local galaxies of similar 
luminosities. Using improved estimates for the IR luminosities of these sources, we investigate whether 
their infrared emission is found to be in excess relative to that expected based on extinction corrected 
UV star formation rates (SFRs), possibly suggesting the presence of an obscured AGN. Through a 
spectral decomposition of mid-infrared spectroscopic data, we are able to isolate the fraction of IR 
luminosity arising from an AGN as opposed to star formation activity. This fraction is only able to 
account for '^30% of the total IR luminosity among the entire sample and ^35% of the "excess" IR 
emission among these sources, on average, suggesting that AGN are not the dominant cause of the 
inferred "mid-infrared excesses" in these systems. Of the sources identified as having mid-infrared 
excesses, half are accounted for by using proper bolometric corrections while half show the presence 
of obscured AGN. This implies sky and space densities for Compton-thick AGN of ~1600 deg~^ and 
~1.3 X 10~^ Mpc"'^, respectively. We also note that IR luminosities derived from SED fitting the 
mid-infrared and 70 /im broadband photometry agree within ^^50% to those values estimated using 
the additional mid-infrared spectroscopic and submillimeter data. An inspection of the FIR-radio 
correlation shows no evidence for evolution over this redshift range. However, we find that the SMGs 
have IR/radio ratios which are a factor of '^3 lower, on average, than what is measured for star-forming 
galaxies in the local Universe. 

Subject headings: galaxies: evolution - infrared: galaxies - radio continuum: galaxies 
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1. INTRODUCTION 

A reliable estimation of the star formation rates 
(SFRs), stellar mass, and active galactic nuclei (AGN) 
activity at all redshifts is required in order to obtain a 
comprehensive understanding of galaxy evolution. Var- 
ious studies have shown that the SFR density increases 
by a factor of '^5—10 between z ^ and z ~ 1 and 
that >60% of the tota l SFR densi ty at z ^ 1 is ob- 
scured by dust (lElbaz et al. 2002; Chary fc Elbaz I I2001I: 
iLe Floc'h et al. I 120051 : iMagnelh et. al. I I2OO90 . Until 
now, these results have relied on deep surveys at 15 fim 
and 24 /xm which sample the redshifted mid-infrared 
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spectral energy distribution (SED) of galaxies. 

The mid-infrared (5 — 40 /im) SED is a com- 
plex interplay of broad emission features thought to 
arise from polycyclic aromatic hydrocarbon (PAH) 
molecules, silicate absorption features at 9.7 and 18 /im, 
and a mid-infrared continuum from very small grains 
(iLeger fc Puget |[l98l lAllamandola. Tielens. k Barker"! 
1 1985( 1. Mid-infrared luminosities measured near ^^8 /tm 
have been found to correlate with the total infrared 
(IR; 8 ^1000) luminosities of galaxi es in the l ocal Uni- 
verse (|Charv fc Elbaz 2001: E lbaz et al. I [2002 ). which 
itself is a measure of a galaxy's SFR. While a correla- 
tion is found, there does exist a large amount of scat- 
ter jPale et al. 1 120051 : ISmith et al.1l2007t lAmius et al. I 
|2007| ) and systematic departures for low metalhcity sys- 
tems ([Engcibracht_et_al. 2005; Mad den et al. Il2006( l. It 
is therefore uncertain if these empirical relations between 
mid-infrared and total IR luminosities (^ir) for galaxies 
in the local Universe are applicable for galaxies at higher 
redshifts. 

Recently, data from deep far-infrared surveys such as 
the Far-infrared Deep Legacy Survey (FIDEL; PI: M. 
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Dickinson) are confirming the same redshift evolution of 
the SFR density seen in the mid-infrared and do not 
show any ev olution in the SEP of infrared luminous 
galaxies (e.g. iMagnelli et. al. |[2009f ). Furthermore, in- 
dividual lensed Lyman break gal axies (LBGs) at z ^ 2.5 
(iSiana et al. 1120081 IWilson et al . 2008: G onza l ez et al. I 
seem to show bolometric corrections similar to 
those of local galaxies. In contrast. i Rigbv et al. I ()2008f ) 
have provided evidence that the high redshift lensed mid- 
infrared selected galaxies might show a factor of ^2 
stronger rest-frame 8 /im emission compared to their 
total IR luminosities. Part of the contradiction might 
be due to the intrinsic scatter in the physical proper- 
ties of the galaxies, particularly the range of dust color 
temperatures and emissivities se en for galaxies in the lo- 
cal Universe (e.g. lDunne et al. 1 12000 ). Individual galax- 
ies detected in the far-infrared or submillimeter at high 
redshift, for which the comparisons are made, might be 
sampling the scatter in the ratio rather than the median 
trend which the empirical galaxy templates follow. 

In order to assess the accuracy of mid-infrared derived 
bolomet ric luminositi e s as a tracer of star-formation at 
2: ~ 2, iDaddi et al. I ()2007al [bl) compared mid- infrared 
(i.e. observed-frame 24 \xvd) based SFRs to those de- 
rived from extinction corrected UV emission and radio 
emission. They found a systematic offset, whereby the 
mid-infrared derived bolometric luminosities were sys- 
tematically higher than alternate SFR estimates, espe- 
cially among the most intrinsically luminous mid-infrared 
sources. Furthermore, they found that stacking the X- 
ray data on these "mid-infrared excess" galaxies revealed 
the presence of a hard X-ray source suggestive of ob- 
scured AGN activity. Since mid-infrared wavelengths are 
sensitive to w arm dust continuum emission arising from 
AGN activity. iDaddi et al.~l (|2007bf) concluded that these 
mid-infrared excesses arise from deeply obscured AGN 
and report a space density of Compton thick AGN which 
is twice that of all X-ray detected AGN at z ^2. 

The detection of a source in hard X-rays does not nec- 
essarily imply that an AGN is dominating its energetics 
as is found for submillimeter galaxies (SMGs) which ap- 
pear to be both star- formation driven and commonly de- 
tected in the X-rays ([Alexander et al. ll2005t[Pope et al. I 
[2008a). Furthermore, evidence for "mid- infrared excess" 
relative to extinction corrected UV luminosities might 
simply be evidence for optically thick star-formation, a 
phenomenon which is ubiquitous among ultraluminous 
infrared g alaxies (ULIRGs; Ltr > 10^^ i©) in the local 
Universe ([Buat et al. II2007I ). Mid- infrared spectroscopy 
of galaxy samples selected by their mid-infrared flux den- 
sities provides a powerful technique to isolate these var- 
ious selection effects among infrared luminous sources. 
By decomposing these spectra, into PAH and continuum 
components, one can get a handle on the relative frac- 
tion of power being emitted by AGN activity within such 
sour ces relative to the output related to st ar formation 
(e.g. ISaiina et al. Il2007t iPope et al. ] '2008a'). 

The Great Origins Observatories D eep Survey-North 
(GOODS-N; [Dickinson et"an l2003bn field, due to the 
wealth of multiwavelength data spanning the ultraviolet 
to radio wavelengths, is the optimal field in which to 
undertake such a relatively unbiased survey. Using mid- 
infrared spectroscopy for 22 sources selected at 24 /^m 
in the GOODS-N field, along with existing 850 /im and 
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Galaxy Positions and Redshifts 
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62 10 47.3 
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62 17 8.1 
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79 
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16 
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62 11 46.2 
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12 37 4.34 
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18 


12 37 11.37 
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1.996 
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19 


12 37 16.59 


62 16 43.2 
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12 37 18.27 
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additional 70 [ivn imagery obtained as part of FIDEL, we 
aim to improve our understanding of the star formation 
and AGN activity within a diverse group of 0.6 z <^ 2.6 
galaxies. This is done through a more proper estimate 
of IR luminosities and the ability to decompose these 
measurements into star-forming and AGN components. 

The paper is organized as follows: In §2 we introduce 
the sample and describe the observations. The determi- 
nation of IR luminosities from SED fitting, calculation of 
SFRs, and the decomposition of AGN and star-formation 
power using the mid-infrared spectroscopy are described 
in §3. The results are presented in §4 while their impli- 
cations are discussed in §5. Finally, in §6, we summarize 
our conclusions. 

2. THE SPECTROSCOPIC SAMPLE AND 
MULTIWAVELENGTH PHOTOMETRY 

The mid-infrared spectroscopic sample presented here 
consists of 22 sources drawn from 24 /xm detections 
in the Great Obse rvatories Origins Deep Survey-North 
(GOODS-N; Dic kinson et al. Il2003bf) field. Galaxy coor- 
dinates and IDs used throughout the text are given in Ta- 
bled! The sample was chosen to span a range in redshift 
from 0.6 z < 2.6 and contain a variety of object types: 
SMGs, X-ray or optically selected AGN, typical lumi- 
nous infrared galaxies (LIRGs; 10^^ < Ljr < 10^^ ^0), 
and optically faint (i.e. zab > 25 mag) 24 /im sources. 
The classification for each object is given in Table [2] To 
ensure high signal-to-noise in the IRS spectra, a flux den- 
sity cut of fi,{2A /im) > 200 ^Jy was made, although the 
median flux density of the sample is 510 /iJy. The mul- 
tiwavelength photometry is given in Table [3l 
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2.1. IRS Spectroscopy 

Spitzer IRS observations were taken as part of Spitzer 
program GO-20456 (PI: R.-R. Chary) in April/May 
2006. These low resohition (i? = A/AA ~ 100) spectro- 
scopic observations were made using the spectral star- 
ing mode, which observed each target at two nod posi- 
tions (per cycle) within the slit. In Table [2] we provide 
the ramp times and the total number of cycles used per 
low-resolution module. We note that a study focusing 
on the SMGs inclu ded in this sample can be found in 
iPope et al. I ()2008af ) where complete details on the ob- 
servation strategy and data reduction steps are also pre- 
sented. The actual spectra, shifted into the rest frame, 
are shown in Figure [1] along with associated errors. The 
expected positions of the 3.3, 6.2, 7.7, 8.6 and 11.3 ^m 
PAH features, along with the the 9.7 /im silicate absorp- 
tion feature, are indicated with vertical lines. 

For cases where PAH features in the IRS spectra are 
clearly ide ntified, redsh if ts were measured directly as de- 
scribed in iPope et al~1 (|2008af) . The redshifts for the 
remaining galaxies were taken from existing optical spec- 
troscopy (Table [T]). We note that there were some dis- 
crepancies in the IRS and optical spectroscopic redshifts 
(i.e. GN-IRS 5,9,15, and 16) in which case the IRS- 
derived redshifts were used. GN-IRS 5 and 9 are both 
SMGs for which the discrepancies between their IRS 
and optical redshif ts are discussed in Appendix B of 
IPope et al. 1 (Hooli). iPope et al. 1 (|2006( l report that the 
optical {z = 0.884) redshift associated with GN-IRS 16, 
also an SMG, is not the true counterpart. The discrep- 
ancy for GN-IRS 15 is thought to arise from a superpo- 
sition of sources. The optical morphology of this source 
is quite complex, having several components and a range 
of colors. Thus, we strongly suspect we are seeing a su- 
perposition of a mid-infrarcd-bright galaxy located at 
z = 2.20 and a foreground galaxy at z = 1.362. 

The final redshift range among the sample spans 0.6 < 
z < 2.6 with a median of 1.88. A total of 15 galaxies have 
a redshift between 1.4 < z < 2.6, which is almost exactly 
the nominal redshift range fo r the BzK co l or selec tion 
where a comparison with the iDaddi et al. I (|2007b[ ) re- 
sults are appropriate. In addition, 7 of the galaxies have 
redshifts between 0.6 ^ z ^ 1.2 and are intended to be 
representative of the LIRG population. 

2.2. Mid- and Far-Infrared Spitzer Imaging 

Infrared imaging of the GOODS-N field was obtained 
at 16, 24, and 70 /im using the Spitzer Space Tele- 
scope. The 24 fim observations of GOODS-N were taken 
as part of the GOODS legacy program (PI: Dickinson) 
and reach an RMS of ~5 ^Jy (see ICharv et al. Il2007f ). 
The spectroscopic sample presented here spans the range 
220 - 1220 fiJy. 

IRS peak-up observations at 16 (|Teplitz et al. I 
[2005h were made covering an area slightly less extended 
than the GOODS-N 24 /im map. Only one of the sam- 
ple galaxies, GN-IRS 1, lies outside the areal coverage 
of the 16 /im map. These 16 /im data have an RMS of 
~6 /iJy and only sources with a signal-to-noise (S/N) ra- 
tio greater than 5 were considered as a detection. Of the 
21 sources in the 16 /im image field-of-view, 2 were not 
detected and their flux densities were set to the upper 
limit value of 30 /iJy. The 16 /im flux densities among 



the detected sources span a factor of '^25, ranging from 
^40 — 1050 /iJy with a median value of ~280 fj,3y. 

MIPS 70 /im observatio ns were carri ed out as part 
of program GO-3325 (PI: iFraver et al~l [2006) and the 
Far Infrared Deep Extragalactic Legacy project (FI- 
DEL; PI: Dickinson). The full field extends beyond 
the 24 /i m coverage and data processing details can be 
found in iFraver et al. I ()2006[ ). The typica l point-source 
noise of the 70 /im data is ~0.55 mJy (jFraver et al. I 
l2006f ). ~1.6 times larger than the confusi on noise level 
of ffc — 0.35 ± 0.15 mJy estimated by IFraver et al. I 
(|2009( l. Sources having a flux density > 2 mJy and 
a S/N > 6 were considered detections, leading to a to- 
tal of 9 sources; all other sources were assigned an upper 
limit of 3 mJy since this value corres ponds to a differen- 
tial completeness of --80% (Magn eUi et. al. Il2009 ). The 
median 70 /im flux density among the detected sources is 
~5.5 mJy and ranges between ~2.2 — 14.2 mJy, spanning 
a factor of roughly 6. The photometric uncertainty in- 
cluding calibration uncertainties and confusion is ~10% 
at 16 and 24 /im and ~20% at 70 /im. 

2.3. X-ray Imaging 

X-ray imaging of GOODS-N was obtained with 
the Chandra X-ray Ob servatory for a 2 Ms exposure 
(| Alexander et al. I (20031 ). The full band (0.5 - 8.0 keV) 
and hard band (2.0 — 8.0 keV) on-axis sensitivities are 
~7.1 X lO"-^'' and 1.4 x 10~-^^erg cm~^ s~^, respectively. 
Of the 22 sample galaxies, 1 2 were matched w ith an 
X-ray counterpart reported bv I Alexander et aT~l ([2003). 
For the remaining 10 sources, upper limits were measured 
assuming P = 1.4. 

Among the 12 X-ray detected sources we find that the 
0.5 — 8.0 keV flux spans more than 2 orders of magni- 
tude, ranging from 9 — 1120 x 10"^^ erg cm~^ s~^, with 
a median flux of 56 x 10^^^ erg cm^^ s^^ (Table [3]). 
Eight sources have secure detections in the hard band 
(2.0 — 8.0 keV) while only upper limits could be mea- 
sured for the remaining 14 sources. Of those sources 
detected, the range in the 2.0 — 8.0 keV flux spans be- 
tween 25 — 1050 X lO"^*" erg cm~^ (i.e. a factor of 
~40) with a median flux of ^100 x 10^^ erg cni^^ s^^. 

2.4. Optical Imaging 

The GOODS-N fleld has been imaged extensively in 
four bands, -B435, Vgoej *775, and zss o using the ACS 
came ra on the Hubble Space Telescope (iGiav alisco et al.1 
120041 ). Additional ground-based imaging of GOODS- 
N in the J7-band was obtained using the prime- focus 
MOSAIC camera on the KPNO Mayall 4 m telescope 
(IGiavalisco et al. II2004I : ICapak et al. Il2004( l. The 3-cr 
upper Hmits of the UBViz photometry are 0.027, 0.048, 
0.048, 0.091, and 0.131 /iJy, respectively. 

2.5. Radio Continuum Imaging 

The GOODS-N fleld was observed with the Very Large 
Array (VLA) for a total of 165 hours using all 4 conflgu- 
rations. (G. Morrison et al. 2008, in preparation). Given 
that the VLA is a centrally condensed array, the observ- 
ing time per conflguration was scaled as follows, A-array 
1), B-array (1/4), C-arr ay (1/16), and D-array (1/64) 
Owen fc Morrison! l2008t ) . This integration time scal- 
ing provides the best sensitivity for larger sources. The 
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TABLE 2 
Spitzer IRS Observations Summary 

Number of Cycles per Module" 





SLl 


LLl 


LL2 


Selection 


ID 


(240 s) 


b (120 s)*" 


(120 s )'' 


Criteria^ 


1 




40 


31 


OF,X-ray,SMG 


2 




10 


6 


SMG 


3 


10 




6 


X-ray, LIRG 


4 




70 




OF,X-ray,SMG 


5 




55 


35 


OF,SMG 


6 




55 




OF.SMG 


7 




45 


35 


OF,SMG 


8 


15 




12 


X-ray, LIRG 


9 




45 




X-ray, SMG 


10 


10 




12 


X-ray 


11 


15 




20 


X-ray,SMG 


12 




40 


31 


OF 


13 






35 


SMG 


14 


12 




8 


X-ray, LIRG 


15 




16 


12 


OF 


16 


20 




12 


OF,X-ray,SMG 


17 




45 




NL-AGN 


18 




28 


20 


X-ray, SMG 


19 




28 


20 


OF,X-ray 


20 




29 


20 


OF 


21 




10 


6 


X-ray 


22 


10 




12 


LIRG 


Note. — " 


Observations 


were made 


using the IRS star- 



ing mode which observes each target at two nod positions 
per cycle. ^ Ramp times per cycle. Selection criteria: 
OF - optically faint; X-ray detected; SM G - submillime- 
ter galaxy (included in lPope et al. I2008al ): LIRG - typical 
LIRG; NL-AGN - narrow-line AGN. 



final radio map has a local RMS at the phase center of 
^4 /iJy. All sources in the IRS spectroscopic sample were 
matched with a S/N> 4 radio source to a radius of 1" 
except for GN-IRS 13, which had a S/N~3.3. The 20 cm 
flux densities span a factor of '^9 (i.e. ~21 — 202 /iJy) 
with a median value of ~98 /iJy. 

2.6. Submillimeter Imaging 

The GOODS-N field was imaged at 850 fim by a 
number of programs using the Submillimeter Common- 
User Bolometer Array (SCUBA) and the data were com- 
bined into the SC UBA "supermap" (jBorvs et al. Il2003t 
iPope et al. 1 120051 ). The 8 SMGs targeted for IRS spec- 
troscopy were selected from the "superma p" sources with 
robust identifications (|Pope et al. Il2006f) in addition to 
3 SCU BA photometry sources from Ichapman et al. I 
dlMl)- Only 2 of the 11 SMGs are detected at 70 ^inl 
(Huynh et al. 2007). The 850 /im flux densties span 
a range between 2 — 9 mJy with a median flux density 
5.2 mJy. For the remaining non-detected IRS targets, 3- 
(T upper limits were calculated from the 850 /im SCUBA 
supermap. The only exception is GN-IRS 20 which lies 
outside the areal coverage of the SCUBA maps. 

3. DATA ANALYSIS 

In the following section we describe our methods for 
calculating infrared (IR; 8 — 1000 /im) luminosities (^ir) 
for each of the sample galaxies along with IR and UV 
(1500 A ) based SFRs. For all calculations we assume a 
Hubble constant of 71 km s^^, and a standard ACDM 
cosmology with flu — 0.27 and J^a — 0.73. Due to the 
small number of sample galaxies, all average properties 



of the sample are given as a median rather than taking 
a mean. 

3.1. Determining Infrared Luminosities from SED 
Fitting 

Infrared luminosities are extrapolated by fitting the 
mid-infrared spectra and photomet ric data to the SED 
templates of ICharv fc Elbazi (|2001h . and then integrat- 
ing between 8 — 1000 ijm. We ch oose to fit to the 
templates of ICharv fc Elbaz I (|2001h since they have 
been found to exhibit 24/70 /xm flux density ratios 
which are more representative of the average actu- 

1 than either the 



iDale & He 


ou 1 ()2002f) or Laeache. 


templates ( 


MaEneUi et. al. 1120091) 



^ The best-fit SEDs are 

determined by a minimization procedure in which the 
SED templates are allowed to scale such that they are 
being fitted for luminosity and temperature separately. 
Consequently, the amplitude and shape of the SEDs scale 
independently to best match the observations. Weights 
are derived using the uncertainties of the input photom- 
etry (i.e. the S/N ratio of each input photometric band). 
When fitting the mid-infrared spectra, having N individ- 
ual measurements for each wavelength in the spectra, the 
individual flux density is weighted by its S/N ratio and 
an additional factor of 1/A^ so as to not give too much 
weight to the mid-infrared portion of the spectrum in the 
fitting. The normalization constant is then determined 
by a weighted sum of observed-to-template flux density 
ratios for all input data used in the fitting. Errors in 
the fits are determined by a standard Monte Carlo ap- 
proach using the photometric uncertainties of the input 
flux densities. 

Since we are interested to see to what extent th e mid- 
infrared excesses reported bv lDaddi et al. I ()2007al f0) can 
be solely attributed to improper SED fitting rather than 
the presence of an embedded AGN, we co mpute IR lu- 
minosit ies in a number of ways. First, like iDaddi et al.~l 
(|2007al ). we per form the SED fit ting using th e 24 /zm flux 
density alone; ICharv fc Elbaz I (|2001h and iDale et al.~l 
(2001) SED templates are fit independently, and their 
integrated 8 — 1000 /im luminosities are then averaged. 
The associated IR luminosity is designated as Lf^. The 
difference between the resultant L^^ from the two sets 
of SED templates is characterized as a systematic error. 
We emphasize here that IR luminosities derived in this 
manner assume that the mid-infrared to total-infrared 
correlations imprinted in local SED templates are valid 
at all redshifts. 

Next, we use all available spectroscopic and photomet- 
ric data in the SED fitting: IRS spectra, 16, 24, 70, and 
850 /im flux densities where possible. IR luminosities for 
4 out of the 22 sources that are neither detected at 70 
or 850 ^im (i.e. GN-IRS 12, 17, 19, and 20) are esti- 
mated by averaging the IR luminosities from the best- 
fit Chary fc Elbaz (2001) and Dale et al. (2001) SEDs. 
Due to the non-detection of the 4 sources at wavelengths 
close to the peak of the far-infrared emission, these IR 
luminosities are considered to be upper limits. In the 
cases where only upper limit values are available for the 
photometric data, they are not incorporated into the 
minimization but are used to reject fits which have flux 
densities above the associated upper limit. The IR lu- 
minosities derived from these fits are labeled as Lf^ and 
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Fig. 1. — The mid-infrared IRS spectra for each of the 22 sample galaxies. Each spectra is shown with associated errors in the rest- frame 
and has been smoothed to the instrument resolution. The expected positions of the 3.3, 6.2, 7.7, 8.6, and 11.3 fim PAH features are 
indicated by dashed vertical lines; the expected position of the 9.7 ^m silicate absorption feature is indicated with a dotted vertical line. 



are assumed to be the best-fit estimate of the true IR 
luminosity for each source. In addition, we also calculate 
IR luminosities by fitting the SED templates to all avail- 
able photometry (i.e. including the 850 fim data where 
possible) to assess how well the PAH features can be 
characterized by fitting photometric data alone. These 
values are labeled as We give Lf^, Lf^, and Lf^°^ 

values for each galaxy in Table |4] 

Since it is more likely that mid-infrared spectroscopy 
and submillimeter imaging will not be available for a 
large number of sources in many deep fields, we also 
calculate IR luminosities as above excluding these data. 
A comparison of these IR luminosities, identified as 
^IR^^'^"' ^'^ those calculated using all available photom- 
etry will help quantify the reliability of SED fitting using 
16 — 70 /im photometry alone. The IR luminosities de- 
rived by this fitting method are also given in Table [4] and 
are considered to be upper limits for sources not detected 
at 70 /xm. 

To illustrate the differences in the best-fitting SEDs 
using these four met hods, we plot the best-fitting 
ICharv fc Elbaz I (|2001l ) templates along with the spec- 
troscopic and photometric data, for 3 different sources 
(GN-IRS 2, 4, and 15) in Figure [H A comparison of 
each fitting method zoomed in on the mid-infrared range 
of 5 — 15 fim, where PAH emission may be present, is 
also shown for each galaxy in Figure [3] along with the ac- 
tual spectra. While the SMGs i n our sample have been 
presented in Pope et al. I ()200 8aV we note that our SED 
fitting method differs from theirs. We do not include an 
extinction term and exclude the radio photometry which 



allows us to assess the FIR-radio correlation among the 
galaxies. C onsequently, for the SMGs in our sample we 
find that the lPope et al. I (|2008af ) luminosities range from 
being a factor of ~0.8 to ~4.6 times that of our derived 
IR luminosities, and are a factor of ~1.2 times larger, on 
average. 

3.2. Infrared Luminosities from the Radio 

In the local Universe there exists a remarkably 
tight correlation between the far-infrared (FIR; 42.5 — 
122.5 /im) and predominantly non-therm al radio contin- 
uum emission from star-forming galaxies (de Jong et al J 

i l98d ^^Hclou, Soifcr, & Rowan-Robinson 198^; 

lYun. Reddv. fc Condon I I200l[ ). Following a similar 
quantitative treatment of the FIR-radio correlation (i.e. 
iHelou. Soifer. fc Rowan-Robinson I Il985( l. we use the 
total IR (8 — 1000 /im) luminosity, rather than the FIR 
fraction, to parameterize the IR-radio correlation such 
that, 

giR = log L „ _ , ) ■ (1) 



3.75 X 1012^^(20 cm) J 

For a samp l e of 1 64 galaxies without signs of AGN ac- 
tivity, iBen] ()2003[ ) report a median qm value of 2.64 and 
a scatter of 0.26 dex. 

We can estimate the rest-frame 1.4 GHz radio lumi- 
nosities of our spectroscopic sample using. 



L^(20 cm) = AttDIS^{20 cm)(l 



(2) 



where _Dl is the luminosity distance. This calculation 
includes a bandwidth compression term of (1 -I- z)~^ and 
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a iiT-correction of (1 + z)" to rest frame 1.4 GHz. This 
assumes a synchrotron power la w of the form Sy (x v~°' 
where the spectral index a ~0.8 (jCondon 1119921 ). 

Thus, by setting < giR 2.64, we can derive the 
infrared luminosities from the radio specific luminosities 
with: 

Lf]f = 1.64 X 10^^^1,(20 cm) (3) 

We assign a factor of ^2 uncertainty to these luminosities 
since this is the intrinsic sc atter among and within (e.g. 
iMurphv et al. I l2006l . 120081 ) star-forming systems in the 
local Universe. The radio based IR luminosities and qir 
ratios are given in Tabic [3] for each galaxy. 

3.3. Estimating Star-formation Rates 
3.3.1. IR Based SFRs 

Using the conversion given in iKennicutt I (|1998f ). cali- 
brated for a ISalpete"r1 (|1955f ) IMF with mass hmits be- 
tween 0.1 and 100 Mq, we can express our IR luminosi- 
ties as SFRs such that 



SFRiR 



yr" 



: 1.73 X 10-1" ( ^ 

V ^0 



(4) 



Since we have estimated IR luminosities in different man- 
ners, we identify each corresponding SFR accordingly; 
SFR24: SEDs fit by 24 ^m flux densities alone; SFR^|: 
SEDs fit using all available data (i.e. our best-fit esti- 
mates); SFR}j^: using the IR-radio correlation to derive 
corresponding IR luminosities for which to estimate a 
SFR. 

3.3.2. UV Based SFRs 

We also estimate SFRs using rest-frame UV spe- 
cific luminosities of each galaxy where possible. Af- 
ter putting the observed UBViz data into the rest- 
frame, 1500 A flux densities were computed by lin- 
early interpolating along the power-law fits between pho- 
tometric data points landing in the wavele ngth range 
between 1250 - 2600 A. A gain takin g the IKennicutt I 
(fl99l) conversion, assuming the same [ Salpeter I (llQSSl ) 
IMF (0.1 - 100 Mq), we can express the resultant UV 
specific luminosities in terms of a SFR such that. 



SFRuv 

Mq yr-1 



= 5.36 X 10^ 



Lp, Hz-i 



(5) 



For a total of 11 galaxies (i.e. half of the sample) there 
are less than two photometric points in the rest-frame 
wavelength range of 1250 - 2600 A. In these cases the 
flux upper limits were used to set a corresponding upper 
limit on each source's rest-frame UV luminosity. 

To estimate the amoun t of extinction at 1500 A we use 
the ni ethod described bv iMeurer. Heckman. fc Calzetti I 
(|1999l ) which relates the amount of extinction to the 
slope of the UV continuum between 1250 — 2600 A , de- 
flned as (3. Since the redshifts of the sample galaxies are 
all less than z ~ 2.6, any correction to the UV slope 
due to absorption by the intergalactic medium is negligi- 
ble. While the relation of lMeurer. Heckman. fc Calzetti I 
(fr999.) is ca librated for Aisoq, we use modeled extinctio n 
curves (i.e. IWeingartner fc Draine"1l2001t IPraine l[2003f ). 
to translate this into an extinction at 1500 A such that 
^1500 ~ l-05Ai6oo leading to a final relation of 



500 



4.65-1-2.09/3. 



(6) 



We designate SFRs estimated using the extinction cor- 
rected UV luminosities as SFRyy. Of the 11 sources 
having less than two photometric detections in the rest- 
frame wavelength range of 1250-2600 A, two (GN IRS 16 
and 20) have at least one detection and one reliable up- 
per limit in this wavelength range allowing us to esti- 
mate a reasonable lower limit for f}. For the remaining 9 
galaxies, a UV extinction correction was not applied (i.e. 
(3 K, —2.2). The observed and extinction corrected UV 
specific luminosities are given, along with the values of 
/?, in Table H 

For the 11 sources which do not have a sufficient num- 
ber of photometric points in the 1250 — 2600 A range to 
properly measure (3, the SFR is considered to be a lower 
limit given that the extinction correction could not be 
measured and is likely significant for these sources; i.e. 
due to the fact that the sample is 24 /Ltm selected, the 
optical non-detections are probably not from the source 
being intrinsically dim, but rather from the extinction 
factor being extremely large. 

3.4. X-ray Luminosities 

To compare IR and X-ray energetics of these systems, 
the rest-frame 0.5 — 8.0 keV and 2.0 — 8.0 X-ray lumi- 
nosities are calculated as 

Lx = ^^DlJ^{l -f zf-^ (7) 

where L'l is the luminosity distance, /x is the X-ray flux 
in a particular band, and F is the observed or assumed 
photon index (Table [3]). The 0.5 - 8.0 and 2.0 - 8.0 keV 
luminosities are given in Table O We also note here that 
3 of the 12 X-ray sources (i.e. 25%) are optically faint, 
having only IRS-derived redshifts between 1.7 < z < 
1.9. This is precisely the redshift ran ge where optically 
faint X-ra y sources are pred icted to lie ([Alexander et al. 1 
,2001; .Mai nieri et 3171120051 ) . 

3.5. Estimating AGN Contributions 

Using the IRS data, we follow t he procedure described 
in detail by iPope et al." * (2008a') to quantify the AGN 
contribution to the total IR energy budget. This is done 
by first decomposing the IRS spectra into two compo- 
nents; one associated with star formation (PAH tem- 
plate) and the other associated with the hot dust con- 
tinuum emission arising from an embedded AGN. The 
star-forming comp onent is fit us i ng th e starburst com- 
posite template of iBrandl et al. I ()2006f ) while the AGN 
component is characterized by a power law with both 
normalization and slope as free parameters. The effect of 
extinction is considered independently for the PAH and 
continuum c om ponen ts usin g a mo deled ex tinction curve 
q.e IWeingartner fc Drair^ 120011 : iDraine M2003I) . Upon 
finding the best fit between the spectra and combination 
of PAH template plus continuum component, the frac- 
tion of mid-infrared emission arising from the continuum 
component is taken to be the AGN contribution. For 
those galaxies whose spectra did not have high enough 
S/N ratios or lacked PAH features to allow for a proper 
decomposition, an upper limit to the AGN strength was 
obtained by fitting a power-law to the continuum of the 
spectra and integrating below it. The mid- infrared AGN 
fractions are given for each galaxy in Table 21 

Next, to characterize the fraction of the IR luminos- 
ity arising from an AGN, we scale the SED template 
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TABLE 4 

Infrared Luminosities and Associated Parameters 



log L24 log L\^^^^^° log LP^°* log L^l log LfjC ^^^^ aGN 

ID {Lq) (Lq) (LqT (Lq5 [Lq) (Le) Fraction (%) giR 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



1 


12.79 


<12.46 


12.12 


12.12 


13.25 


11.89 


0.50 


1.51 


2 


13.62 


12.90 


12.90 


13.04 


13.12 


12.69 


0.79 


2.56 


3 


11.85 


11.79 


11.79 


11.78 


12.00 


< 11.04 


< 0.25 


2.42 


4 


13.49 


<13.00 


12.48 


12.83 


12.79 


12.21 


0.49 


2.69 


5 


12.86 


<12.44 


12.73 


12.73 


13.25 


11.86 


0.46 


2.13 


6 


12.84 


<12.60 


12.73 


12.71 


12.45 


< 11.78 


< 0.48 


2.90 


7 




<12 29 


12 21 


12 21 




11 43 


n 1 s 

U. -LO 


1 62 


8 


11.64 


11.64 


11.64 


11.59 


11.41 


< 10.77 


< 0.23 


2.82 


9 


12.77 


<12.54 


12.09 


12.09 


12.88 


11.63 


0.30 


1.86 


10 


12.03 


11.86 


11.86 


11.89 


11.86 


11.35 


0.44 


2.67 


11 


12.28 


12.39 


12.38 


12.38 


12.81 


< 11.63 


< 0.36 


2.22 


12 


12.93 


<12.48 


<12.48 


<12.33 


12.93 


< 11.83 


< 0.52 


2.05 


13 


11.66 


<11.50 


11.49 


11.50 


11.80 


< 11.34 


< 0.51 


2.35 


14 


11.90 


11.91 


11.91 


11.91 


12.15 


10.99 


0.22 


2.41 


15 


13.55 


12.78 


12.78 


12.76 


12.74 


12.44 


0.69 


2.66 


16 


13.07 


<12.53 


12.39 


12.39 


12.90 


< 11.93 


< 0.51 


2.14 


17 


13.19 


<12.51 


<12.51 


<12.74 


12.56 


12.42 


0.95 


2.82 


18 


13.13 


<12.54 


12.35 


12.35 


13.11 


< 11.80 


< 0.33 


1.89 


19 


12.94 


<12.23 


<12.24 


<12.23 


12.84 


11.69 


0.29 


2.04 


20 


12.97 


<12.47 


<12.47 


<12.29 


12.47 


11.86 


0.41 


2.47 


21 


13.23 


12.65 


12.65 


12.65 


12.87 


12.40 


0.62 


2.43 


22 


11.66 


11.60 


11.60 


11.60 


11.92 


< 11.04 


< 0.29 


2.32 



Note. — Col (1): Source ID. Col (2): 24-fim derived IR luminosity. Col (3): IR luminosity derived from 
SED fitting photometric Spitzer data: 16, 24, and 70 fim photometry where available. Col (4): IR luminosity 
derived from SED fitting all photometric data: 16, 24, 70, and 850 f-LUi photometry where available. Col (5): 
Best-fit IR luminosites derived by SED fitting all available data: IRS spectroscopy, 16, 24, 70, and 850 fim 
photometry where available. Col (6): IR luminosity derived using K-corrected 20 cm flux densities and the 
FIR-radio correlation (see ^T2\ . Col (7): The AGN contribut ion to the total IR luminosity derived using a 
mid-infrared spectral decomposition of the IRS data (see ^3.5| l. Col (8): The mid-infrared ACN fractions (see 
§3.5). Col (9): The logarithmic IR/radio ratio using the usual q-based parameterization given in Equation 

m 



of Mrk 231 (a nearby ULIRG wh ose IR luminosity is 
known to be dominated by an AGN; lArmus et al. |[2007l ) 
to the AGN fraction of the rest-frame '^5 — 12 fj.m lumi- 
nosity (i.e. the luminosity in the rest-frame wavelength 
range having IRS coverage per source) and integrate the 
SED from 8 - 1000 /im. This is illustrated in Fig ure |4] 
where we show the best-fit ICharv fc Elbaz I (|200lh SED 
for GN-IRS 7 along with the observed IRS spectra, the 
scaled SED of Mrk 231 (i.e. the AGN contribution), and 
the difference between the best-fit SED and the scaled 
SED of Mrk 231 (i.e. the contribution by star forma- 
tion alone). The observed mid-infrared spectra clearly 
matches the shape of the best- fit SED. 

The values of Lfjf ^ are given for each galaxy in Ta- 
ble m By subtracting these values from the IR lumi- 
nosites obtained through our SED fitting, we can esti- 
mate the amount of IR luminosity arising from star for- 
mation alone such that 



/ SFRfg^GN 
V Mq yr-i 



= 1.73 X 10" 



L 



all,noAGN 
IR 



^0 



(8) 



where Lf°^'^^ 



rail 



7- AGN 
'^IR ■ 



Similarly, we can also cal- 
culate an AGN-subtracted SFR from the 24 /Ltm data by 
SED fitting the 24 /im photometry after first subtracting 
off AGN fractions extrapolated from those given in Table 
m the corresponding IR luminosities and SFRs are de- 

, J r24,noAGN j ot^t-, 24.noAGN i ^t,t 

noted as L^^ and SIRjj^ , respectively. We 

have assumed that the underlying mid-infrared contin- 
uum arises purely from an AGN, however hot dust asso- 
ciated with vigorous star formation may also contribute 
to the mid-infrared continuum emission as well. There- 



fore, the mid-infrared AGN fractions, and associated es- 
timates for the IR luminosities of the AGN, should all be 
considered conservative upper limits. 

4. RESULTS 

We now present the results for the various SED fitting 
methods among the sample while looking for trends with 
redshift. In doing so, we also compare the infrared-based 
SFRs to those estimated from extinction corrected UV 
emission for galaxies having redshifts between 1-4 < z < 
2.6. 

4.1. Accuracy of Derived IR Luminosities 

In F igure [5] we plot the best-fitting ICharv &: Elbaz I 
(|200l[) SEDs for 3 sources (GN-IRS 2, 4, and 15) to il- 
lustrate the three fitting methods. While the addition of 
the IRS spectra and submillimeter data to the available 
Spitzer 16, 24, and 70 fim photometry does not seem 
to significantly alter the choice of the fit, the fits using 
24 /im flux densities alone are highly discrepant. Zoom- 
ing in on the 5 — 15 nm mid-inf rared regions of the four 
different ICharv &: Elbaz I (|2001h fits for each galaxy in 
Figure [31 one easily sees that by using 24 fim fiux densi- 
ties alone the equivalent widths of the aromatic features 
observed in the IRS s pectra are underestimated by the 
ICharv fc Elbaz I ()2001[ ) templates, resulting in highly dis- 
crepant estimates of IR luminosity. 

The range of IR luminosities derived by fitting the 
24 /im photometry alone span a range between 4.3 x 
10^^ — 4.2 X IG^'^ ^0 (i.e. nearly two orders of magnitude) 
with a median of 7.5 X 10^^ Lq. Looking at the best-fit es- 
timates for the IR luminosities among the entire sample. 
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TABLE 5 
X-RAY AND UV Luminosities 



log io.5-8.0 kcV log i2. 0-8.0 kcV log ^UV log iyV 

ID (Lq) (Lq) (Lq Hz-1) (Lq Hz-1) 13 

(1) (2) (3) (4) (5) 



1 


10.16 


10.13 


<-5.50 


<-5.50 




2 


< 8.91 


< 9.14 


-4.74 


-3.23 


-0.47 


3 


8.14 


< 8.03 


<-4.94 


<-4.94 




4 


9.59 


9.49 


-5.11 


-3.09 


0.14 


5 


< 8.75 


< 8.94 


<-5.43 


<-5.43 




g 


< 8.85 


< 9.07 


<-5.32 


<-5.32 




7 


< 8.64 


< 8.84 


<-5.44 


<-5.44 




8 


7.96 


7.86 


<-4.78 


<-4.78 




9 


8.67 


8.72 


-4.85 


-3.46 


-0.61 


10 


9.79 


9.48 


<-4.41 


<-4.41 




11 


8.69 


< 8.47 


-4.83 


-2.96 


-0.03 


12 


< 8.81 


< 9.04 


-4.98 


-3.76 


-0.82 


13 


< 7.85 


< 8.05 


-5.20 


-3.18 


0.14 


14 


8.12 


< 8.16 


-6.18 


-3.69 


0.71 


15 


< 8.98 


< 9.24 


-5.18 


-4.04 


-0.91 


16 


8.43 


< 8.56 


<-5.59 


<-4.04 


<-0.42 


17 


< 8.80 


< 9.05 


-4.23 


-3.96 


-1.95 


18 


8.97 


8.94 


-4.57 


-2.53 


0.18 


19 


8.75 


8.73 


<-5.53 


<-5.53 




20 


< 9.24 


< 9.49 


<-5.50 


<-2.69 


< 1.09 


21 


9.23 


9.20 


-5.55 


-2.48 


1.40 


22 


< 7.87 


< 8.07 


<-5.60 


<-5.60 





Note. — Col (1): Source ID. Col (2): K-corrcctcd full band (0.5 - 8.0 koV) lumi- 
nosities. Col (3): K-corrcctcd hard band (2.0 — S.OkeV) luminosities. Col (4): Linearly 
interpolated 1500 A specific luminosities. Col(5): Extinction corrected 1500 A specific 
luminosities. Col(6): Slope of t he UV continuum between 1250 — 2600 A which is used 
to estimate extinction foUowing lMeurer. Heckman. Calzetti I ( I1999D . 



we find that they range between 3.2 x 10^^ — 1.1 x 10^^ Lq 
(i.e. a factor of ~35) having a median of 2.2 x 10^^ Lq. 
The 24 /xm-derived IR luminosity is a factor of ~4 larger 
than that for the best-fit estimates, on average. 

Using our best-fit IR luminosites, we also note that 16 
out of the 22 sample galaxies are classified as ULIRGs 
while the remaining 6 are LIRGs. Looking at the differ- 
ences between the 24 /xm to our best-fit IR luminosities, 
the discrepancy is found to be much larger among the 
ULIRGs than for the LIRGs. The 24 /.tm-derived IR lu- 
minosities are larger than the best-fit IR luminosities by 
a factor of on average, for the ULIRGs while being 
only ~1.2 times larger, on average, for the LIRGs. 

4.2. Strength of PAH Features 

To quantify how well the mid-infrared spectral features 
are matched by the various SED fitting methods, we com- 
pare the ratio from the peak of the 7.7 /im PAH feature 
to the median flux density in the 9 — 11 /zm trough re- 
gion (i.e. 9.7 /zm silicate absorption feature). Given that 
the wavelength coverage varies for each object, making 
it difficult to properly es timate the true continuum levels 
(e.g. iPope et al. Il2008a[) , we choose to compare this ratio 
rather than actual PAH equivalent widths. Furthermore, 
choosing these features enables us to calculate this ratio 
for nearly the entire sample; we do not have spectro- 
scopic coverage over the rest-frame 9 — 11 fj,m region for 
GN-IRS 16. 

When fitting the SEDs using the 24 iim photometry 
alone the peak-to-trough ratios are ~0.24±0.15 times the 
peak-to-trough ratios measured for the observed spectra 
(i.e. a factor of ^4 times smaller), on average. By in- 
stead fitting the SEDs using all of the photometric data, 
the peak-to-trough ratios are now ~0.75±0.84 times the 
ratios of the actual data (i.e. ^35% smaller), on aver- 



age. Similarly, when the mid-infrared spectra is included 
in the fitting with all of the photometric data, the peak- 
to-trough ratios are found to be ~0.75 ±0.58 times the 
ratios for the the observed spectra, on average, suggest- 
ing that fitting with mid- and far-infrared photometry 
alone does a reasonable job choosing templates with the 
correct mid-infrared spectral features. Consequently, it 
appears that the improper estimates of IR luminosities 
from fitting with 24 fim photometry alone are most likely 
due to the fact that local high luminosity galaxies do not 
exhibit large PAH equivalent widths relative to galaxies 
of similar luminosity a t higher reds h ifts. This result is 
consistent with that of iDesai et al. I (|2007| ) who showed 
that local ULIRGs typically have PAH equivalent widths 
that are significantly smaller than those measured for 
ULIRGs at high-redshift. Since the strongest PAH fea- 
tures, lying in between 6-12 /im, begin to redshift into 
the observed 24 /im band at z > 1.4, we expect the dis- 
crepancies between the fitting methods and estimates for 
the IR luminosity of galaxies, to increase with increasing 
redshift. For the sample galaxies in the redshift range 
between 1.4 ^ z < 2.6, we find that the peak-to-trough 
ratios of the SEDs chosen by fitting the 24 /im data alone 
and using all of the available photometric data are an av- 
erage factor of ~0.19 ± 0.11 (i.e. ~5 times smaller) and 
^^0.88 ± 0.94 (i.e. ~14% smaller) times those measured 
for the data, respectively. 

We also note that by remeasuring the peak-to-trough 
ratios of the best-fit SEDs after subtracting the fitted 
AGN component, the median ratios are ^1.11 ± 0.30 
and ^1.49 ± 0.42 for the entire sample and galaxies at 
z > 1.4, respectively. This result demonstrates that 
the AGN/star-formation decomposition is fairly reliable 
and that, as previously stated, we likely overestimate the 
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true AGN contribution by assuming that all of the mid- 
infrared continuum emission arises from hot dust associ- 
ated with the AGN. These statistics exclude 3 galaxies 
(GN-IRS 1, 13, 17) for which the AGN component over- 
subtracts the continuum of the best-fit SEDs resulting in 
a negative peak-to-trough ratio. Including them into the 
above calculations does not affect the median values, but 
increases the dispersion significantly. 

4.3. Redshift Dependency 

In Figure [H we plot Lf^, Lf^, and versus redshift 
for the 70 /.tm detected galaxies. Although the different 
IR luminosity estimates largely agree for the z < lA 
sources, the 3 sources at z > 1.4 have IR luminosities 
derived from fitting the 24 /im flux densities alone that 
are larger by an average factor of ^^4 than what is found 
when fitting the SEDs with all available data. 

Since the 70 /xm detections could potentially be sam- 
pling objects with the warmest FIR dust color tempera- 
tures, we also compare the different IR luminosity esti- 
mates of all sources in our sample in Figure [6l We find a 
significant change in this ratio with redshift. While the 
objects near z ~ 1 all have ratios around unity, galaxies 
in a redshift range between 1.4 ^ z < 2.6 exhibit ratios of 
~4.6 ± 1.4, on average. We also note that galaxies lying 
at redshifts below and above z ~ 1.4 typically have 
values below and above ~3 x 10^ ^ Lp,, respectivel y . This 
result is consistent with that of iPapovich et al~l ()2007f ) 
whom report that IR luminosities derived from stacking 
70 and 160 data for a sample of 24 /xm-bright (i.e. 
/i,(24 ^m) > 250 ^Jy) galaxies lying in a redshift range 
between 1.5 z < 2.5, are a factor of 2 — 10 lower than 
those inferred from using only 24 iim flux densities. A 
similar comparison with the radio continuum based lu- 
minosities over the same redshift range shows that the 
24 /im based IR luminosities are higher by a factor of 
~1.5 ± 1.9, on average. 

It is unclear if this discrepancy is due to the flux limit 
of our sample which implies that we are sampling sources 
with high luminosities at z ~ 2 for which no local analogs 
exist. This would suggest that the extrapolation of the 
local SEDs to these high redshift sources is invalid. A sec- 
ond possibility is that z ^ 2 sources have sufficiently low 
metallicities that the infrared SEDs deviate from those 
of z ~ 1 or local galaxies. A third possibility that has 
been suggested is that warm dust continuum emission 
from an obscured AGN in these sources is biasing the 
24 /i m fiux upward and re sulting in overestimates of the 
SFR (|Daddi et al. ll2007aD . 

4.4. AGN-Subtracted SFR Estimates from 24 fJ-ni 

The SFRs inferred from the 24 /im-derived IR lumi- 
nosities overestimate those inferred from our best-fit IR 
luminosities, corrected for any contributions from AGN 
(SFRf^^^^; see §3.5) by a factor of ~6, on average, 
among the entire sample. To see whether this overes- 
timation by fitting the 24 /im photometry is the result of 
excess mid-infrared emission arising from hot-dust emis- 
sion associated with the presence of an AGN, we refit 
the 24 /im fiux densities to the SED templates after cor- 
recting them by the corresponding mid-infrared AGN 
fraction. We find that SFRs calculated using IR lumi- 
nosities derived by SED fitting 24 /im flux densities af- 
ter subtracting the mid-infrared AGN contributions are 



^2.5 ± 2.3 times larger, on average, than the values of 
SFRfp^'^'^ among the entire sample. There also appears 
to be a redshift dependence such that the SFRs calcu- 
lated using IR luminosities derived by SED fltting 24 /im 
flux densities after subtracting the mid-infrared AGN 
contributions are '-^1.0 ± 0.4 and ~3.3 ± 2.2 times larger, 
on average, than the values of SFRJ"^^*^^ for galaxies at 
redshifts below and above z ~ 1.4, respectively. 

We believe that IR luminosities are overestimated by 
fltting local SEDs to 24 /im photometry alone because of 
differences in PAH equivalent widths among galaxies of 
similar luminosity at different redshifts. This conclusion 
is supported by our findings in §4.2 where we compared 
the observed PAH strengths to those from the different 
SED fits. Alternatively, this result could also imply that 
using Mrk 231 to approximate the AGN component is 
inappropriate. However, given that we believe our mid- 
infrared AGN fractions and IR luminosity estimates for 
the AGN to be conservative upper limits, this is unlikely 
to be a dominant effect. 

4.5. Multiwavelength Comparison between SFR 
Estimates 

In order to assess the reliability of the different SFR 
indicators, we make a cross comparison between the UV, 
radio and IR SFR estimates. We examine how these 
estimates compare for galaxies located at redshifts below 
and above z ~ 1.4. 

4.5.1. z < 1.4 Galaxies 

For the 7 galaxies located at z < 1.4, we find that SFRs 
inferred from their 24 /im data alone range from being 
^0.8 to ~1.4 times that of those SFRs calculated using 
our best-fit IR luminosities, and are only 16% larger, on 
average. However, using the radio data to infer a SFR via 
the FIR-radio correlation, we find that the radio SFRs 
overestimate those calculated using the best-fit IR lumi- 
nosities by a factor of ~1.7, on average, and range from 
being a factor of '^0.7 to ~2.7 times that of the SFRs 
measured using are best-fit IR luminosities. Among these 
7 galaxies, only 3 have measurable UV slopes allowing 
us to calculate a UV corrected SFR. These SFRs range 
from being a factor of ~0.8 to ^^6.5 times that of the 
SFRs measured using our best-fit IR luminosities. 

4.5.2. z > 1.4 Galaxies 

iDaddi et al. I ()2007a[ ) have identified discrepancies be- 
tween the UV, radio, and, IR derived SFR estimates in 
z ~ 2 galaxies selected using the BzK selection tech- 
nique. We select the 15 galaxies in our sample with 
1-4 z < 2.6 which matches the redshift range spanned 
by the BzK selection. Figure [7| shows the the ratio of 
SFRs inferred from the 24 /tm data (plus the uncorrected 
UV-based SFR) to those estimated from extinction cor- 
rected UV measurements for these 15 galaxies. Galaxies 
which are not detected in the optical, such that their UV 
slopes could not be estimated, are plotted as upper lim- 
its although they are difficult to interpret. Their best-fit 
IR-derived SFRs are nearly 2 orders of magnitude larger 
than the uncorrected UV-based SFRs, on average. How- 
ever, their extinction corrected UV SFRs could poten- 
tially be very large, reducing the plotted ratio. 
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iDaddi et all (|2007af) defined sources for which 

to be "mid-infrared excess" sources. While the numer- 
ator adds the contribution of energetic emission from 
young massive stars that is not absorbed and reemitted 
by dust grains to obtain a total SFR, the extincted UV 
SFR is often negligible compared to the IR-based SFR 
term. 

Not surprisingly, we find that all 7 galaxies for which 
we do not have UV slope measurements, and thus can- 
not make UV extinction corrections, all meet the mid- 
infrared excess criterion. For these sources to fall below 
the mid-infrared excess criterion would require an aver- 
age extinction of Ay ^1.7 mag; and to fall to a ratio near 
unity requires an average extinction of Ay > 2.2 mag. 
As for the remaining 8 sources, we find 6 are identified 
as having mid-infrared excesses: GN-IRS 2, 4, 9, 12, 15, 
adn 17. 

The two sources for which we do not find any mid- 
infrared excess, due to their steep UV corrections, are 
GN-IRS 18 and 21. Both GN-IRS 18 and 21 are X-ray 
detected, however GN-IRS 18 is also an SMG. Morpho- 
logically, GN-IRS 18 appears extended possibly due to a 
major merger event. The extinction estimated for GN- 
IRS 18 is Ay ~ 1.92 mag, only 0.16 mag larger than the 
median extinction found for all sample galaxies for which 
the rest-frame UV slope could be measured. The extinc- 
tion estimated for GN 21, however, is the largest among 
the entire sample at Ay « 2.89 mag. While their opti- 
cal extinctions differ significantly, these are the only two 
sources having extinction corrected UV SFRs in excess 
of 1000 Mq yr-i. 

By plotting again the ratio of IR -I- UV to UV corrected 
SFRs, this time instead using our best-fit IR luminosi- 
ties, versus redshift in the middle panel of Figure [3 we 
find that 2 of the 6 previously identified mid-infrared 
excess galaxies (GN-IRS 4 and 9) are no longer mid- 
infrared excess sources. These two galaxies are also the 
only ones among the 6 to have hard X-ray detections. In 
fact, only 2 sources, GN-IRS 2 and 15, seem to remain 
definitive "mid-infrared excess" galaxies given that GN- 
IRS 12, which barely lies above the mid-infrared excess 
threshold, and GN-IRS 17 are not detected at either 70 
or 850 /zm implying that their associated best-fit IR lu- 
minosities are technically upper limits. If, on the other 
hand, we use the radio continuum based SFRs to con- 
struct the same plot, as shown in the bottom panel of Fig- 
ure [3 we find that while points have shifted downward, 
only GN-IRS 4 is no longer considered a mid-infrared ex- 
cess source leaving 5 of the original 6 mid-infrared excess 
galaxies above the threshold. 

To summarize, we find that by using the best-fit IR 
luminosities we can definitively account for 2/6 mid- 
infrared excess sources (i.e. ^ of those found using 24 /xm 
derived IR luminosities). This fraction is increased to | 
by assuming that we are still overestimating the IR lumi- 
nosities of the two galaxies for which we only have upper 
limits at far-infrared wavelengths. The median ratio of 
best-fit IR to extinction corrected UV SFRs for sources 
in the redshift range between 1.4 < z < 2.6 is ~4. By 
instead using the radio continuum and FIR-radio cor- 
relation to derive IR luminosities, the median ratio be- 
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Fig. 2.— The best-fit SEDs of lCharv fc Elbaz I 1 I2OOII) for GN- 
IRS 2 (top panel), 4 (middle panel), and 15 (bottom panel) using 
only 24 ^lra flux densities {dotted lines), all 3 Spitzer (16, 24, and 
70 fj,m) flux densities {dashed lines), all Spitzer and 850 fim pho- 
tometry {dot- dashed lines), and all Spitzer photometry along with 
the IRS spectra and the 850 fim flux densities {solid lines). The 

best-fit SEDs associated with Lj^'^'*'™, and Lf]^ for GN- 

IRS 15 are very similar, which is why the dashed and dot-dashed 
lines overlay each other. 

tween the radio derived-to-extinction corrected UV SFRs 
is even larger at ~7. We investigate in the discussion be- 
low whether the contribution from an obscured AGN is 
large enough to account for the discrepancy between the 
IR and UV SFR estimates among these sources. 

5. DISCUSSION 

Using IRS spectroscopy along with additional photo- 
metric data, we have set out to characterize better the IR 
properties for a diverse sample of 24 /xm selected sources. 
In doing this we also try to explain the true nature of 
sources considered to have an excess of IR emission rel- 
ative to their ultraviolet luminosity when 24 fim pho- 
tometry, alone, is used to estimate SFRs. Specifically, 
through a spectral decomposition of the IRS data, we 
determine that half of the " mid- infrared excess" sources 
are Compton-thick AGN. A summary of our findings for 
each of the 22 sample galaxies is given in Table [6l 

5.1. Discrepancies in SFR Estimates versus the 
Contribution from A GN 

As discussed in i) l4.5.2i we find that even after using the 
best-fit IR luminosities, Lf^, to derive SFRs, we still find 
4/8 sources in the redshift range between 1.4 < z < 2.6 
with measurable UV slopes to exhibit excesses in the 
infrared relative to what is measured by extinction cor- 
rected UV data (middle panel of Figure [7|) . To see if 
the observed excess is in fact due to the presence of an 
obscured AGN, we estimate the fraction of AGN power 
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Fig. 3.— A blowup of the 5 - 15 range for the best-fit SEDs of Chary & Elbaz (200l!) for all sources using only 24 fim flux 
densities {dotted lines), all 3 Spitzer (16, 24, and 70 fim) flux densities {dashed lines), Spitzer and submillimeter (850 fira) photometry 
(dot-dashed lines), and all photometry {Spitzer and submillimeter) along with the IRS spectra {solid lines). The actual spectra, smoothed 
to the instrumental resolution, are overplotted in red. The SEDs chosen when fitting with the 24 fim data alone rarely characterize the 
observed PAH emission as compared to when longer wavelength data is included in the fitting. 
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Fig. 4.— The best-fit SEDs of ICharv fc Elbaz I l(200ll ) for GN- 
IRS 7 plotted as the bold solid line, given in the rest-frame. Over- 
plotted are the IRS spectra of the source {solid histogram- style 
line), the scaled Mrk 231 template {dashed line), corresponding to 
the contribution to the SED by the AGN, and the contribution to 
the SED by star formation {dot-dashed line). 



contributing to the total IR luminosity using our mid- 
infrared spectroscopy as described in ij3.5l 

In the left panel of Figure [8] we plot the ratio of our 
best-fit IR derived SFRs, plus an unobscured UV com- 
ponent, to the extinction corrected UV SFRs against the 
mid-infrared AGN fractions before and after subtracting 
the AGN contribution to the total IR luminosity. We find 
that 3/8 galaxies having redshifts between 1.4 S z < 2.6 
and which continue to exhibit an excess of IR emission 
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Fig. 5. — We plot IR luminosity, estimated by three different 
means, versus redshift for galaxies in our IRS sample which have 
70 /xm detections; associated 1-cr error bars are given. Ljj| val- 
ues {open circles) are calculated by SED template fitting using the 
24 fim flux density only while values {filled circles) also in- 
clude the 16, 70, and 850 pm photometry (where available) along 
with IRS spectra. values are calculated using the well-known 

FIR-radio correlation and plotted as open squares; since the uncer- 
tainties associated with these estimates are large (i.e. a factor of 
~2), we do not plot error bars. The Lj^ estimated by SED fitting 
the sample's /iy(24 fim) flux density limit of 220 fijj as a function 
of redshift is plotted as the dashed line. 



have mid-infrared AGN fractions that are larger than 
60%. Using the Mrk 231 SED as outlined in § 3.5, we 
can estimate the contribution of AGN to the total IR 
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Fig. 6. — We plot the ratio between the 24 /xm-derived IR 
luminosity (L^^) to that derived using the additional constraints 
from the IRS spectroscopy along with the 16, 70, and 850 fim 
data. Sources having hard band (2.0 — 8.0 keV) X-ray detections 
are identified with an "X". For sources not identified at 70 and 
850 /im, we plot an upward arrow indicating that the IR luminosity 
is only an upper limit. 

luminosity. 

Looking at the entire sample (i.e. all 22 galaxies), we 
find that AGN account for ~30% of the total IR lumi- 
nosity, on average, with a range of '^10 — 70% (Table|4]). 
This is similar to the average AGN contribution of ~30% 
we find among the hard band X-ray detected sources. In- 
terestingly, we also find that the average contribution to 
the total IR luminosity by AGN is also '^30% for both 
sub-samples of SMGs and non-SMGs. This contribution 
by the AGN to the IR luminosity of SMG s is a factor 
of ^^2 times larger than that reported by iPope et al. I 
(|2008af) which we attribute to the differences in the IR 
luminosities arising from differ ences in the SEP fitting 
methods; as previously stated, iPope et al. I (|2008 a) in- 
cluded the radio photometry in their SED fitting leading 
to IR luminosities typically larger than those presented 
here. 

In the right panel of Figure [8] we plot AGN luminosity 
versus the difference between the 24 /im and best-fit es- 
timate IR luminosities. A clear trend of increasing AGN 
luminosity with increasing overestimation of the IR lu- 
minosity using 24 /im photometry alone is observed; or, 
in other words, AGN luminosity appears to scale with 
increasing mid-infrared luminosity. Performing an ordi- 
nary least squares fit to the data, we find that 

/ r AGN \ / r 24 _ r all 

log i^^j = (0.50±0.04) log ' i« 



-(5.61±0.52), 
(10) 



While the differences between the 24 /im and best-fit 
IR luminosities arc large for galaxies having large (i.e. 
>60%) mid-infrared AGN fractions, we find these quan- 
tities to be unrelated for galaxies having smaller mid- 
infrared AGN fractions. We also find that the AGN 
luminosity accounts for only 16% of the difference be- 
tween the 24 fim derived and best- fit IR luminosities, 
on average. This indicates that the AGN contribution 
to the mid-infrared excesses, which would result in over- 
estimates of IR-based SFRs, is close to negligible when 
compared to the improper bolometric correction applied 
when estimating total IR luminosities by fitting local 
SED templates with 24 /im data alone. As shown in § 
4.4, SFRs derived by subtracting the AGN fraction from 



the 24 fim photometry are larger than the AGN corrected 
best-fit IR SFRs by a factor of ^^1 to ~3.3, on average for 
redshifts below and above z ^ 1.4, respectively. Further- 
more, the AGN-corrected 24 /im SFRs are factors of ~1 
to ^68 times larger than the corresponding UV-corrected 
SFRs for sample galaxies at redshifts above z ~ 1.4 . 
Therefore, the AGN is not the dominant source of the 
inferred "mid-infrared excesses" among these systems. 

Since the FIR-radio correlation is thought to solely 
arise from the processes associated with massive star for- 
mation, we looked to see if the IR/radio ratios themselves 
are sensitive to the fractional AGN power by comparing 
the ratio of /Lf^ with giR values. A trend was not 
found suggesting that giR is insensitive to the fractional 
IR output by the AGN. This result likely arises because 
the AGN contributes to the radio and infrared emission 
by different amounts. 

5.1.1. The Persistence of "Mid-Infrared Excess" Sources 

By subtracting off our estimate of the AGN contribu- 
tion to the total IR luminosity, we determine whether 
this is enough to account for the remaining IR excesses 
for galaxies in Figure \7\ We recreate the top and mid- 
dle panels of Figure [7] in the top and bottom panels of 
Figure [51 respectively, by computing the 24 /im-derived 
and best-fit IR-based SFRs after subtracting off the esti- 
mated AGN contribution for the 8 galaxies having optical 
detections which allowed for proper estimates of extinc- 
tion corrected UV SFRs; the decrease in IR SFRs are 
indicated by vertical lines. 
Among these sources we find that (Lj^ — 



^IR 



1/L 



IR, 



are 



both ~45%, on average. However, even by using the 
AGN-corrected 24 /im SFRs, none of the previously 
identified mid-infrared excess sources are removed. 
While the AGN fractions are similar in both cases, the 
AGN-subtracted 24 /im SFRs are still ~6 times larger, 
on average, than the AGN-corrected SFRs using the 
best-fit IR luminosities. In the bottom panel of Figure 
[9] we find that by subtracting the AGN luminosities we 
are able to account for the observed IR excesses in only 
one more source; GN-IRS 12, which is plotted as an 
upper limit, is now no longer considered to have excess 
IR emission. Therefore, of the 8 sources with optical 
detections such that we can estimate a UV extinction 
correction in the redshift range between 1.4 ^ z < 2.6, 6 
of which would be classified as "mid-infrared excesses" 
using 24 /im-based IR SFRs, 3 (GN-IRS 2, 15 and 
17) remain mid- infrared excess sources even after 
subtraction of the AGN component. 

Ultimately, we find that even after correcting for the 
presence of an AGN, the IR-based SFRs are a factor of 
~2.8 times larger, on average, than those derived from 
extinction corrected UV measurements for 1.4 < z < 2.6 
galaxies with secure estimates of their rest-frame UV 
slopes. Or, in other words, AGN are only able to ac- 
count for ^35% of the excess IR emission measured, on 
average. Since we do not believe our IR luminosities are 
overestimated by such a large amount, given the small 
uncertainties associated with our SED fitting, we inter- 
pret this result to suggest that using the UV slope to 
derive an extinction correction is not appropriate for the 
entire sample; extinctions are underestimated thereby 
yielding underestimates for the UV-derived SFRs. This 
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Fig. 7. — In the top panel we plot the ratio of SFRs estimated from the 24 /im flux densities (plus a non-obscured UV-derived SFR) 
to extinction corrected UV SFRs as a function of redshift. The dashed line indicates where the logarithm of this rat io is equal to 0.5 dex; 
galaxies taking values larger than this are considered to be "mid-infrared excess" sources (i.e. IDaddi et af~ll2007alV Sources undetected 
in the optical, such that its UV slope could not be fit, are plotted with downward arrows. In the middle panel we plot the same ratio 
versus redshift, except that we use aU available data to derive the IR based SFRs (SFRf^). The IR based SFRs are considered to be upper 
limits for those galaxies not detected at 70 and 850 ^m The same is done again in the bottom panel, except that we use the 20 cm radio 
continuum flux densities and the FIR-radio correlation to derive SFRs (SFR||^). Mid-infrared excess sources largely remain even when the 
bolometric corrections are accounted for. 



may not be too surprising given that there are catas- 
trophic failures using the rest-frame UV slope to derive 
the extinction for galaxies having both low and high ex- 
tinctions. If extinction within a galaxy is too high, this 
method will fail simply because the galaxy's ISM will be- 
come optically thick. In the case of low extinction the 
UV slope will no longer be as sensitive to the amount of 
extinction since variations arising from contributions by 
a galaxy's old stella r population will become important 
(jReddv et al. l[2006l ). 

5.1.2. AGN Dominated Sources 

Of the remaining 3 mid-infrared excess sources after 
correcting for the contribution of the AGN to the total 
IR luminosity, GN-IRS 17 is the only source plotted as an 
upper limit, however all have mid-infrared AGN fractions 
>60% suggesting that they are AGN dominated. While 
the remaining IR excess for GN-IRS 17 may arise from 
improper SED fitting due to only having upper limits at 
70 and 850 /im, we note that its optical morphology is 
point-like and consistent with being an AGN. It may then 
be that the SED templates used for the fitting, which are 
based on star- forming galaxies, are not appropriate re- 
sulting in an incorrect estimate of the total IR luminos- 
ity. In fact, this source has recently been i dentifie d as a 
Compton-thick AGN bv [Alexander et al. I (| 2008aV Sim- 
ilarly, GN-IRS 2 is t hought to be a Compton-thick AGN 
(|Pope et al. Il2008al ). so it is not surprising to find that it 



lies above the threshold even though it is not plotted as 
an upper limit. Thus, the only remaining mid- infrared 
excess source which is not plotted as an upper limit or 
has not been previously identified as a Compton-thick 
AGN is GN-IRS 15. 

Inspecting the IRS spectra for this source, it shows 
strong evidence of a deep silicate absorption feature and 
weak PAH equivalent widths, suggestive of a deeply em- 
bedded AGN. The ACS imagery of this source is also 
point-like, which too seems consistent with the source 
being an AGN. Consequently, we believe that GN-IRS 15 
is primarily powered by an AGN (i.e. a Compton-thick 
AGN) for which our SFR estimates are not physically 
meaningful. 

The other source which appears to be clearly AGN 
dominated according to its mid-infrared AGN fraction is 
GN-IRS 21. We note that this source never appears to be 
a mid-infrared excess source even when comparing its UV 
corrected SFR to that derived from SED fitting 24 /im 
photometry. However, by attempting to correct for the 
AGN in this source, we find that the best-fit IR derived 
SFR becomes ^5 times smaller than that of the UV cor- 
rected SFR; clearly the UV and IR SFR estimates for this 
AGN dominated source are not physically meaningful. 
We also note that GN-IRS 4 and 12, which just barely 
drop below the mid-infrared excess threshold in the bot- 
tom panel of Figure [9l have mid-infrared AGN fractions 
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Fig. 8. — In the loft panel we plot the ratio of IR-derived SFRs using our best-fit IR luminosities, plus an unobscured UV component, 
to extinction corrected UV SFRs versus the mid-infrared AGN fraction before and after subtracting the AGN contribution to the total IR 
luminosities using filled circles and open diamonds, respectively. The dashed line indicates where the logarithm of this ratio is equal to 
0.5 dex; galaxies taking values larger than this are considered to be "mid-infrared excess" sources (i.e. iDaddi et al. iL2007ai) . In the right 
panel we plot the estimated AGN luminosity (see i|3.5l l as a function of the difference between the 24 fim and best-fit IR luminosities. 
Sources for which only upper limits to the AGN mid-infrared luminosity fraction could be determined are shown as downward arrows. 
An ordinary least squares fit is overplotted as a solid line while the dotted line indicates a one-to-one relationship to show that the AGN 
luminosity is not able to account for the difference in the bolometric correction when fitting the SED templates with 24 fim photometry 
alone. 



near 50%. However, since our mid-infrared AGN frac- 
tions are considered to be upper limits, we believe these 
systems to likely be dominated by star formation. 

Of the 6 originally identified "mid-infrared excess" 
sources, based on extrapolating IR-based SFRs from 
SED fitting 24 /im photometry alone, 3 remain even 
after accounting for the presence of AGN and apply- 
ing proper bolometric corrections. This suggests that 
the sky a nd space densi t ies of C ompton-thick AGN re- 
ported by iDaddi et al. I (|2007bf ) likely overestimate the 
true values by a factor of ~2, yielding corrected values 
of '-1600 deg-2 and ~1.3 x lO""* Mpc'^, respectively 

5.1.3. Are "Mid-Infrared Excess" Sources DOGs? 

In Figure ini we label sources which would be identified 
as dust-obscu r ed ga laxies (DOGs) using the definition 
of lDev et al.l ([2001 . /^(24 ^m)//^(i?) > 982, to see if 
they clearly separate out as mid-infrared excess sources. 
We estimate R band flux densities by interpolating be- 
tween the V and 1775 band photometry. More than half 
of the sample (i.e. 13 sources, 7 of which we could not 
measure proper extinction corrected UV SFRs and are 
therefore not included in Figure [9]) fit the DOG defini- 
tion. Not all DOGS can be classified as mid-infrared 
excess sources. Conversely, not all mid-infrared excess 
sources are DOGs. This result is not too surprising since 
DOGs are thought to be powered by both star formation 
and AGN processes (jPope et al. 1120080 ). However, we 
do find that DOGs are only identified for sources wtih 
z > 1.4. In fact, 87% of our sample having a redshift 
greater than 1.4 are classified as DOGs. 

5.2. Application For Other Fields with Deep Spitzer 
and Herschel Far-Infrared Imaging 

We compare the results of our SED fitting to situations 
when only mid-infrared photometry (e.g. at 16 and/or 
24 /im) along with additional measurements at 70 //m 
are available. Such a comparison is especially important 
for the FIDEL fields, because of their deep 70 /im imag- 
ing, as well as for upcoming dedicated i/ersc/ieZ programs 
such as GOOm-Herschel (GOODS-H; PI: D. Elbaz). 



In the third panel of Figure [10] we plot the ratio of 
IR luminosities derived by fitting the 16, 24, and 70 /im 
photometry to Lfp as a function of redshift. The median 
ratio has a value of ^^1.02 ± 0.51. We find the ratio to 
be approximately constant with increasing redshift, with 
average values of ^1.01 ± 0.06 and '^1.20 ± 0.61 for red- 
shifts below and above 2: ~ 1.4, respectively, suggesting 
that the systematic variations introduced by the pres- 
ence of strong PAH emission redshifing into the 24 //m 
band at z > 1.4 is largely taken into account by having 
the additional 70 ^m data point. We also see that the 
dispersion significantly increases with redshift. We there- 
fore believe that using mid-infrared photometry with the 
addition of a 70 /im measurement provides fairly reliable 
estimates for the true IR luminosity compared to when 
fitting with the 24 /zm data alone. We also note that it 
is the addition of the 70 /.tm data, and not the additional 
16 /.tm photometry, which is responsible for the improved 
determination of IR luminosities. 

In the first (top) panel of Figure [10] we plot the ratio 
of IR luminosities derived by fitting the 16 and 24 /xm 
photometry to those using all of the available data ver- 
sus redshift. The median ratio of 16 -I- 24 /im-derived 
to our best-fit IR luminosities is '^1.17±1.65. As with 
the comparison made above, the median and dispersion 
in the ratio of IR luminosity estimates increases with 
redshift, going from ~0.99±0.19 to -1.29 ± 1.82 at red- 
shifts below and larger than z ~ 1.4. This is of course 
not surprising since fitting the 16 and 24 /im photom- 
etry for redshifts beyond z ^ 1 is complicated by the 
presence of PAH and silicate absorption features being 
redshifted into both bands. Instead, taking the ratio of 
24 -I- 70 /im-derived to our best-fit IR luminosities, we 
find a median ratio of ~1.07±0.45 which is illustrated 
in the second panel of Figure 1101 Similarly, we find the 
median ratio to remain near unity for galaxies below and 
above 2; ~ 1.4, however the increase in the dispersion be- 
tween these subsets is less significant being ^1.04 ±0.13 
and '^1.29±0.53, respectively. Furthermore, we find that 
the additon of the 16 /im photometry to SED fitting the 
24 and 70 /im data does not result in largely different 
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Fig. 9. — The top panel is similar to that of the top panel of Figure [T] except that we have subtracted the AGN fraction from the 
24 fim photometry before calculating the associated IR luminosities and SFRs. We also only show galaxies for which a proper extinction 
corrected UV SFR could be calculated. The decrease, which is ~45% on average among these sources, is shown by a vertical line. In the 
bottom panel we subtract the AGN contribution to our best estimates of the IR luminosity (i.e. Li^ values estimated from SED fitting all 
available data to derive SFRs). The decrease in SFR by correcting for the AGN (~45% on average among these sources as well) is again 
shown by a vertical line. While the decrease to the SFRs by correcting for AGN is similar among these sources, the AGN-subtracted 24 /im 
SFRs are still '^6 times larger, on average, than the AGN-corrected SFRs from our best-fit IR luminosities. 



IR luminosity estimates; the median ratio of IR lumi- 
nosities from fitting the 16 — 70 /im data to those from 
fitting only the 24 and 70 fim data is ^1.03 ± 0.29 and 
ranges from ~1.03 ± 0.12 to '--^l.OB ± 0.34 for galaxies be- 
low and above z ~ 1.4, respectively. Such a comparison 
is useful since most deep fields only have coverage at 24 
and 70 /im. Consequently, the addition of the 16 /im 
photometry doesn't appear to improve the fitting signif- 
icantly compared with the addition of the 70 ^m data. 

It is interesting to see how including the mid-infrared 
spectra affects the SED template fitting compared to 
when using only the available photometry. We plot the 
ratio of Lf^°'/Lfl as a function of redshift in the bot- 
tom panel of Figure [TOl We find that the ratio has a 
median that is approximately unity and a dispersion of 
±0.22. While the median of the ratio remains near unity 
for redshifts below and above z ~ 1.4, we again find that 
the dispersion increases with redshift, taking values of 
~0.06 and ~0.26 in these redshifts bins, respectively. 

Furthermore, we find that using the 70 /itm data alone 
to derive IR luminosities does not appear to be reliable as 
this method significantly underestimate the true IR lu- 
minosity. Among the 70 fim detected sources (only 2/9 
which are SMGs), the best-fit IR luminosities are larger 
by an average factor of ~1.9±0.22 than the IR luminosi- 
ties derived from the 70 fim data alone. Such a result 
is expected if the average galaxy at high redshift has an 
SED which is cooler than that of a a local galaxy at the 



same luminosity (e.g. P ope et al. 1 120061 : iKovacs et al. I 
2006; Huvnh et al. 200'^ However, to make a conclu- 
sive statement about such an effect requires a precise 
sampling of the full IR SED, including the peak, which 
we do not have. 

Seeing that the addition of 70 fim data, alone, to mid- 
infrared photometry yields relatively reliable estimates of 
IR luminosity among these galaxies, FIDEL observations 
(i.e. deep 70 /im imaging) of the entire GOODS-N field 
are used to better constrain estimates of the IR luminos- 
ity for all 24 fim detected sources having spectroscopic 
redshifts in E.J. Murphy et al. (2009, in preparation). 
By doing so, we are able to improve estimates for how 
the IR luminosity density evolves with redshift and char- 
acterize better the number density of mid-infrared excess 
sources. 

5.3. The FIR-Radio Correlation 

Using these data also enables us to look at the behavior 
of the FIR-radio correlation among a diverse group of 
galaxies spanning a redshift range between 0.6 ^ z < 
2.6. In calculating the radio based IR luminosities and 
associated SFRs we have assumed that the FIR-radio 
correlation remains constant with redshift. We will now 
try to verify if such an assumption was appropriate. 

5.3.1. Lack of Evolution with Redshift 
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TABLE 6 

Summary of Sources 



ID 


Class" 


Comment 


^ 


AGN 


KJ V Oi'JJL't; iliU-CLCl lliilidLrt; 


2 


CT ACN 


IVlltl mil dl CU. t;A.ljt;00 CVCli diLrd xlVJlN O U U LI CLL, LiUli dlUJ. UlwUtJi UtJiUlliL Li IL, t^Ui i tjL, LiUii 


3 


SF 


O Li Uiltl i Ji.J.i J_ t,arLLlit;o* KJ V 01LJL/U> llivill- LCJ. illlliCLLC. 


4 


SP 
or 


Strong PAH fG3,tu.rGS j U V-slopG u.iidGrGstirn.3jti6s Gxtrinction 


o 


SP 
or 


UV-slopc iiidctGriniiiSytrG 


a 
o 


SP 

O-T 


UV-slopG iiidctGrminSytrG 


7 


O-T 


Strong PAH fG3,turGS5 U V-slopG indGtGrmin&tG 


g 


SF 


O Li Uiig iTxlli itJdL Lli t;0, KJ V OiUUC iiiLlC LCi liiiiid LC 


g 


SF 


TT\/— Glr»i~»f5 rT\W(^r''i'l'\r PQ+imp+fG PV+inffioTi 

\J V OlUL/t; L-WiiCV_.Liy l_-oLiiila»LCD CALlll^ LiUll 


10 


SF 


T T \/ — g1 r»i~»P 1 TiiH fitr'T'TTi 1 n a tti 

VJ V OlUL/t; IliVi.C Lll-J. illlli(XLC 


11 


SF 


Strong PAH features^ UV-slopG correctly estimates extinction 


12 


SF 


Strong PAH features; UV-slope underestimates extinction; low qi^i value 


13 


SF 


UV-slope overestimates extinction 


14 


SF 


UV-slope correctly estimates extinction 


15 


CT-AGN 


Mid-infrared excess even after AGN subtraction and proper bolometric correction 


16 


SF 


Strong PAH features 


17 


CT-AGN 


Mid-infrared excess even after AGN subtraction and proper bolometric correction 


18 


SF 


Strong PAH features; UV-slope overestimates extinction 


19 


SF 


Strong PAH features 


20 


SF 


UV-slope indeterminate 


21 


AGN 


UV-slope Sz, 24 fim conspire to produce consistent SFRs 


22 


SF 


Strong PAH features; UV-slope indeterminate 



Note. — " Final classification: AGN - IR AGN fraction larger tlian 50%; CT-AGN - Compton-thick AGN; 
SF- star formation dominated source 



In Figure [TT] we plot qir, defined in Equation [TJ as 
a function of redshift. Infrared luminosities are esti- 
mated from SED fitting using all available photometry 
and the IRS spectroscopy. The median qi^ value, plotted 
as dot-dashed line in Figure [TTl is 2.41 ± 0.39 dex, below, 
but c onsistent with, the local value of 2.64 ± 0.26 ( Bell 1 
l2003f ). Interestingly, we note that the 4 galaxies having 
the largest mid-infrared AGN fractions (GN-IRS 2, 15, 
17, and 21) have giR, values very near the canonical ra- 
tio, being within 0.18 dex, on average. This indicates 
that the FIR-radio correlation is a poor discriminant of 
obscured AGN activity. 

While the correlation appears to remain linear with 
redshift for the few objects being studied here, more than 
50% of the SMGs in this sample (i.e. 7/11) have IR/radio 
ratios < 2.24. This is a factor of >2.5 below the aver- 
age value measured in the local Universe for both nor- 
mal star- forming galaxies and ULIRGs. These galaxies 
drive the large dispersion found among the entire sam- 
ple relative to what is found in the local Universe. The 
SMGs show a scatter in giR of 0.44 dex compared to 
0.27 dex for the non-SMGs suggesting a more hetero- 
geneous population. In contra st, 12/13 ULIRGs in the 
lYun. Reddv. &: Condon I (|2001l ) sample, (i.e. excluding a 
single, radio loud AGN), have a median IR/radio ratio 
which is 33% larger than the canonical value and a dis- 
persion of only 0.28 dex. While the average qir value for 
the non-SMGs is 2.43 dex, nearer to the local value of 
2.64 dex, the median qir, value among all SMGs in the 
sample is 2.14 dex, a factor of ~3.2 times lower than the 
canonical value. 

The errors in Lf^ appear too small to explain this 
systematic departure from the canonical IR/radio ratio 
among the SMGs; by increasing Lf^ values by the 3(7 er- 
rors and recomputing giR, the median value among the 
SMGs is 2.16 dex, still a factor for ~3. times lower than 
the canoni cal ratio. Thi s resul t is consistent with the 
findings of iKovacs et al. I (|2006( l who reported similarly 



low IR/radio ratios for a sample of «15 SMGs out to 
z ~ 3.5 detected at 350 and 850 /xm. The possible rea- 
sons for the observed systematic departure off the FIR- 
radio correlation for this galaxy population are discussed 
below. 

5.4. Physical Explanations for Differences in qi^ 
among the SMGs 

Assuming that the uncertainty in our estimates of IR 
luminosity are not responsible for the low IR/radio ratios 
among the SMGs, there are a number of physical expla- 
nations for why galaxies may have IR/radio ratios that 
are significantly lower than the local correlation value. 
First we note that a single spectral index was used to in- 
correct all of the observed radio flux densities. While 
we have assumed a radio spectral index of a « 0.8, 
it would take a negative slope (i.e. a ~ —0.4) to re- 
turn the average IR/radio ratio back to the canonical 
value for the SMGs. Flat or inverted spectra associated 
with such indices are indicative of AGN (e.g. gigahertz 
peaked sources arising from synchrotron self absorption). 
The radio spectral slopes for z ~ 2 ULI RGs are typ- 
ically not found to exhibit such indices; ISajina et al. I 
(|2008| ) report only 3/48 Spitzer selected ULIRGs have 
radio spectral indices between 610 MHz and 1.4 GHz (i.e. 
o^i 4'ghz^) < 0.4. Each of these sources also show very 
weak PAH emission which is unlike the SMGs presented 
here. 

iKovacs et al. I (|2006f ) attributed the low IR/radio ra- 
tios to a possible decrease in the effective dust emis- 
sivity index for SMGs compared to local star-forming 
galaxies. However, this is hard to prove without a fine 
sampling of the full IR SED. The simplest explana- 
tion is extra radio emission due to the presence of an 
AGN. Optical spectroscopy of the SMGs indicate the 
presence of AGN acti vity (e.g. ISwinbank et al. 1 120041 : 
IChapman et al. Il2005l ). Analysis of the X-ray data on 
SMGs also indicates the prese nce of an AGN in most 
SMGs (| Alexander et al. 1120051 ). Thus, it is very likely 
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Fig. 12. — The estimated AGN contribution to the total IR 
luminosites are plotted against the 0.5 — 8.0 keV luminosities. Up- 
per limits limits to the AGN and X-ray luminosities are identified 
by appropriately directioned arrows. The solid line indicates the 
canonical relation between the IR and X-ray luminosities of AGN- 
dominated sources (i.e. quasars) reported by Elvis ct al. (1994). 
The dashed line indicates the fit to AGN-classificd SMGs presented 
in Alexander et al. (2005) and assumes an AGN-IR luminosity 
fraction of 100%. The 3 Compton-thick AGN (i.e. GN-IRS 2, 15, 
and 17) have dust covering fractions significantly larger than that 
of typical quasars. 



Fig. 10. — In the top panel we plot the ratio between the 

16 24 TO 

16 and 24 /jm-derived IR luminosity (ijj^' ' ) to that derived 
using the additional constraints from the IRS spectroscopy and 
850 fira submillimeter data. Similarly, in the second and third 
panels, wc plot the ratio of the 24 and 70 /im-derived (ijj^'^") 

and the 16, 24, and 70 /im-derived (ij^''^^'™) IR luminosities to 
those estimated from fitting all available data, respectively. In the 
bottom panel we plot the ratio between IR luminosities derived 
from fitting all available photometric data (ijf^°' to our best-fit 
IR luminsoties, which also included the IRS spectra in the SED 
template fitting. This comparison is shown to see how well the mid- 
infrared spectral region is characterized by fitting photometry alone 
and how differences are reflected in the associated IR luminosity 
values. Errors in the fitting are shown. In each panel a solid line 
is plotted indicating a ratio of unity. 
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Fig. 11. — We plot q^i ratios, as defined in Equation [T] cal- 
culated using IR luminosities from SED fitting all available data 
against rcdshift. The solid line indicates the local average qi^ value 
of 2.64 (Bell 2003) and the dotted lines are the ztlo" values. The 
thick dot-dash line indicates the median value for the sample and 
falls just above the — 1-cr line. 



that some of the radio emission arises from the nuclear 
black hole as can be seen by the fact that a number of 
the most deviant qm values are in fact detected in the 
hard band (2.0 - 8.0 keV) X-rays. 

In order to assess if the departure from the canonical 
<7iR value for the SMGs is due to radio emission from an 
AGN, we recompute qm ratios by first subtracting the 
AGN component of the IR luminosity and attempting to 
make a similar correction to the radio. This is done by 
assuming that, like the IR emission, the radio emission 
of Mrk 231 is dominated by its AGN. This is reason- 
able since Mrk 231 has a = 2.39, nearly a factor of 
~2 below the nominal value for star-forming galaxies. 
Taking the 20 cm flux density of M rk 231 of 0.309 Jy 
(jCondon. Cotton, fc Broderick Il2002f ) along with its red- 
shift (z = 0.042), the corresponding 20 cm specific lumi- 
nosity is 3.36 x 1O~"^L0 Hz~^. The fraction of radio power 
associated with the AGN for each galaxy was then esti- 
mated by scaling the radio luminosity of Mrk 231 by the 
mid-infared AGN fractions of the mid-infrared luminos- 
ity obtained in H3.5I This contribution is then subtracted 
from the measured 20 cm specific luminosities for each 
galaxy before computing the new IR/radio ratios. 

Rather than decreasing the dispersion among the sam- 
ple galaxies, this attempt to remove emission arising from 
the AGN has actually increased the dispersion from 0.49 
to 0.54 dex. This result suggests that the fractional in- 
frared output of the AGN does not scale with the frac- 
tional output in the radio, consistent with what was 
found in §5.1 by comparing the fractional AGN power 
with giR. 

Another explanation for the low IR/radio ratios 
may be related to galaxy-galaxy interactions. The 
morphology of SMGs seem to suggest major mergers 
which are driving inte nse bursts of star formation (e.g. 
ICh^ man et al~l l2003l ). In the local Universe, a num- 
ber of interacting galaxy pairs (i.e. so called "Taffy" 
galaxies) exhibit IR/radio ratios which are a factor 
of ~2 lower than the canonical value (jCondon et al I 
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[Tool iCondon. Helou. fc Jarrett l[200l . These systems 
are characterized by a synchrotron "bridge", a region 
of bright radio continuum connecting the galaxy pairs 
which is thought to arise from a recent interaction and 
Hkely contain both cosmic rays and magnetic fields. The 
excess radio emission arising from the gaps between these 
galaxy pairs typically accounts for half of the total ra- 
dio emission from the entire system. If SMGs have in- 
deed undergone a similar major-merger event, leading to 
molecular gas and synchrotron bridges between galaxy 
pairs which do not yet facilitate a significant amount of 
ongoing star formation, this scenario may explains the 
lower IR/radio ratios. 

On the other hand, perhaps the extreme star-bursts oc- 
curring within SMGs lead to physical conditions which 
affect the acceleration and cooling processes of cosmic- 
ray electrons which boost the synchrotron emissivity. 
The sample of SMGs presented here is far too small to 
properly address this question, which will be discussed 
in a future paper. 

5.5. IR and X-ray Luminosities 

A comparison of the AGN and 0.5 — 8.0 keV lumi- 
nosities are shown in Figure [T^] for the entire sample. 
There exists a loose relationship between the IR and 
full band X-ray luminosities of AGN dominated sources 
(i.e. quasars) suc h that Lq.s-s.o kov/-^iR ~ 0.025 (e.g. 
lElvis et al. 111994 ) as indicated by the solid line in Fig- 
ure [121 We also include the fit to AG N-classified SMGs 
presented in [Alexander et alTI ([20.0 S) having an average 
io. 5-8.0 kcv/iiR ratio of -0.002. 

When the AGN luminosities are plotted against the 
full band X-ray luminosities we find that nearly every 
galaxy, save GN-IRS 1 and GN-IRS 10, remain either 
near or above the AGN-classified SMG relation and well 
away from the canonical AGN relation line for quasars. 
Despite the fact that GN-IRS 1 and GN-IRS 10 are 
hard band X-ray detected, they have mid-infrared AGN 
fractions which are <50%. In fact, the median mid- 
infrared AGN fraction among all hard band X-ray de- 
tected sources is ~44% with a range of ~23 — 62%. This 
indicates that a hard band X-ray detection does not nec- 
essarily imply a galaxy is AGN dominated. 

We believe that the reason most sources lie generally 
above the above the AGN-classificd SMG relation is pri- 
marily from the overestimation of AGN luminosities; as 
previously stated, the mid-infrared AGN fractions and 
associated AGN luminosities should be thought of as up- 
per limits since hot dust contributing to the mid-infrared 
continuum emission may be powered by star formation. 
Alternatively, there is some uncertainty on the absorp- 
tion corrections to the X-ray flux, and it is possible that 
they c ould be a factor of ~2 higher ([Alexander et al. I 
[2008b( l. It is also possible that the AGN in these sources 
contribute significantly more to the tot al IR output than 
the few percent suggested for SMGs bv [Alexander et al. I 
([2005) arising from a large AGN dust covering factor. 
The 4 galaxies having the largest mid-infrared AGN frac- 
tions (GN-IRS 2, 15, 17, and 21), three of which are 
thought to be Compton-thick AGN (GN-IRS 2, 15, and 
17), all lie well away from the typical quasar relation sug- 
gesting that they have significantly larger dust covering 
fractions than quasars. 



6. CONCLUSIONS 

In the present study we have used observations from 
the mid-infrared to the submillimeter to properly char- 
acterize the IR luminosities for a diverse sample of 22 
galaxies spanning a redshift range of0.6<z<2.6. In 
addition, we have used the mid-infrared spectra of these 
sources to estimate the fractions of their IR luminosi- 
ties which arise from an AGN. Our conclusions can be 
summarized as follows: 

1. IR (8 - 1000 ^m) luminosities derived by SED 
fitting observed 24 /im flux densities alone are 
well matched to those when additional mid-infrared 
spectroscopy and 16, 70, and 850 /im photometry 
are included in the fits for galaxies having z < 1.4 
and Lj^ values typically <3 x 10^^ Lq. In contrast, 
for galaxies lying in a redshift range between 1.4 
and 2.6 with Lf^ values typically >3 x 10^^ Lq, IR 
luminosities derived by SED template fitting using 
observed 24 /xm flux densities alone overestimate 
the true IR luminosity by a factor of ~5, on aver- 
age, compared to fitting all available data. A com- 
parison between the observed mid-infrared spectra 
with that of the SEDs chosen from fitting 24 /im 
photometry alone and from fitting all available pho- 
tometric data demonstrates that local high lumi- 
nosity SED templates show weaker PAH emission 
by an average factor of ~5 in this redshift range 
and do not properly characterize the contribution 
from PAH emission. 

2. After decomposing the IR luminosity into star 
forming and AGN components, we find the AGN 
luminosity to be increasing with increasing differ- 
ence between the 24 /im-derived and our best-fit 
IR luminosities. Such a trend suggests that the 
AGN power increases with mid-infrared luminos- 
ity. However, we also find that the median frac- 
tion of the AGN to the difference between the 
24 /im-derived and best-fit IR luminosities is only 
16% suggesting the AGN power is almost negligible 
compared to the bolometric correction necessary to 
properly calibrate the 24 /im-derived IR luminosi- 
ties. 

3. Infrared luminosities calculated using both mid- 
infrared and 70 /im photometry agree to within 
^^50%, on average, to those calculated using addi- 
tional mid-infrared spectroscopic and submillime- 
ter data and do not show any systematic devia- 
tions with increasing redshift. Including the sub- 
millimeter photometry in this calculation improves 
the agreement between the two sets of IR luminosi- 
ties to within ^^20%. This result implies that IR 
luminosities derived for deep fields having 70 /im 
data, like the FIDEL fields, should be able to ob- 
tain fairly reliable estimates for the total IR lumi- 
nosity of sources up to z ~ 3 even in the absence 
of submillimeter data. This will be investigated for 
GOODS-N in a forthcoming paper. 

4. Even after correcting for the presence of an AGN, 
IR-based SFRs are still larger than those derived 
from extinction corrected UV measurements by a 
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factor of ~2.8, on average, for the 1.4 z < 2.6 
galaxies with secure estimates of their rest-frame 
UV slopes. Or, in other words, the AGN emission 
is only able to account for ~35% of the "excess" IR 
emission, on average, and is not the dominant cause 
of the mid-infrared excesses observed in these sys- 
tems. This suggests that either the SED fitting 
is overestimating the total IR luminosity, which 
seems unlikely given the small uncertainties asso- 
ciated with our fitting, or perhaps that the UV 
extinction correction underestimates the true ex- 
tinction factors for these galaxies. 

5. By using proper bolo metric corrections and cor- 
recting for the presence of AGN, we are able to 
account for half of the sources which are identi- 
fied to_JiavB_a^' rmd4iifrared excess", as defined 
by iDaddi et aTl (|2007ar ). based on their 24 /xm- 
derived IR luminosities. This indicates that the 
sky and sp ace densities of Comp ton-thick AGN 
reported by iDaddi et aTl (|2007b( ) are likely high 
by a factor of ^2. We therefore report cor- 
rected sky density for Compton-thick AGN of 
~1600 deg~^, and a corresponding space density 
of -1.3 X lO-'i Mpc-3. 

6. We do not see any clear signatures of evolution 
in the FIR-radio correlation with redshift out to 



z ~ 2.6 for this sample of galaxies. However, the 
median IR/radio ratio measured for the SMGs in- 
cluded in our IRS selected sub-sample is 2.14 dex, 
a factor of ~3 lower than what is found locally for 
star- forming galaxies. 

7. A comparison of the estimated AGN and full band 
(0.5 — 8.0 keV) X-ray luminosities indicate a higher 
IR output than observed canonically for quasars for 
the majority of X-ray detected sources in the sam- 
ple. This can be explained if the AGN luminosities 
are being overestimated due to star formation con- 
tributing significantly to the hot dust continuum in 
the mid-infrared or if the AGN dust covering factor 
is large in these sources. 
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TABLE 3 

MULTIWAVELENGTH PHOTOMETRY 





r /in — 17 
/(). 5-8.0 kcV/J-U 
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/2. 0-8.0 kcv/10 












/,.(16 fim) 


/^(24 fim) 


/^(70 fim) 


/,y(850 fim) 


8^(20 cm) 
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(erg cm s ) 


(erg cm s ) 


pa 


ipJy) 


ipJy) 


(^Jy) 


(^Jy) 


(^Jy) 


(^Jy) 


(M-ly) 




(/^Jy) 


1 


1120.0 


1050.0 


0.44 


< 0.027 


< 0.048 


< 0.048 


< 0.13 




374 


< 3000 


5400 


201.89 


2 


< 19.9 


< 34.1 


1.40t 


0.032 


0.212 


0.353 


0.71 


496.8 


1220 


4030 


7900 


127.46 


3 


41.0 


< 32.1 


1.45 t 


2.572 


3.993 


11.941 


35.84 


774.7 


1210 


11100 


< 24200 


153.07 


4 


102.0 


79.7 


0.97 


< 0.027 


< 0.048 


0.072 


0.19 


60.2 


303 


< 3000 


4900 


34.33 


5 


< 14.0 


< 21.4 


1.40t 


< 0.027 


< 0.048 


< 0.048 


< 0.13 


42.1 


330 


< 3000 


7500 


171.46 


6 


< 14.4 


< 23.9 


1.40t 


< 0.027 


< 0.048 


< 0.048 


< 0.13 


< 30.0 


215 


< 3000 


5200 


21.72 
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< 10.9 


< 17.4 


1.40+ 


< 0.027 


< 0.048 


< 0.048 


< 0.13 
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347 


< 3000 
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32.1 


25.4 
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3.661 
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113.3 
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< 3000 
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95.91 


10 
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13.29 
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42.69 


11 


28.9 
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1.63 * 


0.752 
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6.33 
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1.40t 
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< 3000 
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Note. — '^Photon index: ' indicates the assumed index used when only upper limits were measured; * indicates that the given photon index is an upper or lower Umit. 



